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Abstract. Anomalous diffusion is predicted for Brownian particles in inhomogeneous 
viscosity landscapes by means of scaling arguments, which are substantiated through 
numerical simulations. Analytical solutions of the related Fokker-Planck equation in 
limiting cases confirm our results. For an ensemble of particles starting at a spatial 
minimum (maximum) of the viscous damping we find subdiffusive (supcrdiffusive) 
motion. Superdiffusion occurs also for a monotonically varying viscosity profile. We 
suggest different substances for related experimental investigations. 



1. Introduction 

Brownian particle dynamics plays a key role for many transport processes in various 
disciplines. Even though a century has passed since the publication of Einstein's seminal 
work on normal diffusion [JJ, the field of diffusion processes still attracts enormous 
interest [21 El HI [5] . The essential signature of normal diffusion is the linear temporal 
growth of the mean square displacement of Brownian particles, 



with a — 1, whereas for a ^ 1 the notion anomalous diffusion has been introduced 
[HI [H El [10]. Being not only of fundamental but also of great technological relevance, 
the dynamical processes underlying anomalous diffusion phenomena have started to 
be investigated during recent decades in systems as diverse as collective ordering 
phenomena |Tll [T2] . particle diffusion in mesoscopic systems [131 03] > or social networks 
|15j . At the same time, progress in single particle manipulation and detection on the /im 
and sub-/im scale fostered investigations on random particle motions [161 E3 HH1 fT9| [20] . 

Diffusion processes with an exponent < a < 1 in (JXJ) are called subdiffusive. 
Several examples of such dynamic behavior are known for systems with a short 
inhomogeneity length scale, such as in amorphous semiconductors [10], groundwater 
motion (21], diffusion in gels [IE], and in several biological systems such as in the 
bacterial cytoplasm [IT], during protein diffusion through cell membranes [18], or as 
a possible method to detect the microstructures in actin filament networks [J_9] . Besides 
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subdiffusive motion, heterogeneous parameter landscapes are considered to be important 
also in generalized reaction-diffusion systems and for bifurcations in pattern forming 
systems in general p2J [22] [23j El]. Conversely, superdiffusive processes are random 
motions with an exponent a > 1, which are found in rather different areas of natural 
and social sciences. Superdiffusion tends to occur in active and driven systems where 
random steps are interrupted by intermediate and nearly deterministic motion, such as 
for particles in (turbulent) random velocity fields [25], [26] , during cell migration [20], for 
active transport in cells [13], or during the spread of diseases [T5j . 

Another class of Brownian motion that attracted remarkable interest recently is 
described by a nonlinear Langevin equation with a velocity dependent viscous damping 
coefficient which is used as a model for active and relativistic motion [271 1211 12H1 [30] as 
well as in the context of ratchet models [31] . 

Here, we investigate the effects of inhomogeneous viscosities on the one-dimensional 
dynamics of Brownian particles, considering viscosity variations, which are slow on the 
scale of the particles' random steps. In this case, the interplay of a diffusing particle with 
its inhomogeneous surrounding leads to intermediate regimes of anomalous diffusion, 
similar to the intermediate Rouse regime of diffusion of polymer segments [32]. Near a 
viscosity minimum we find subdiffusion, whereas superdiffusive behavior is found near 
a viscosity maximum and for monotonously varying viscosity profiles. In all three cases 
the particles' probability distribution is found to be non-Gaussian. 

Our analysis is based on a nonlinear Langevin equation and the appendant Fokker- 
Planck equation, which are presented in section [2] together with the spatially varying 
model viscosities considered in this work. In section [3] we calculate scaling formulas for 
the mean square displacement and analytical solutions of the Fokker-Planck equation 
for several limiting cases. These results are compared with numerical simulations 
in section HJ Finally, we comment on possible experimental realizations within our 
summary in section [5j 



2. Model 

2.1. Langevin and Fokker-Planck equation 

The one-dimensional Langevin equation for the velocity v (t) of a particle of mass m 
and radius a immersed in a medium with a spatially varying viscosity landscape rj(x) is 
given by 

v(t) = -j(x)v + g(x)£(t) , (2) 

with j(x) = 67Tar](x)/'m as the associated damping. The right hand side includes white 
noise characterized by 

(£(*)> = and m at')) = S(t-t'). (3) 

For the noise strength g(x) we assume a local fluctuation-dissipation relation [33] 

= 2 7 (s)fc B T 
y ' m 



Anomalous diffusion in viscosity landscapes 



3 



with fog being the Boltzmann constant and T denoting the constant temperature. In 
the overdamped limit, i. e. for times longer than the characteristic time 7 _1 , a Langevin 
equation for the particle's position can be derived from (j2J) by using the method of 
adiabatic elimination [3j El |3l]. In Ito's interpretation it takes the form 

k B T d 1 



x(t) 



m dx j(x 




(5) 



The corresponding Fokker-Planck equation for the probability density of a particle 
P(x,t) can be deduced straightforwardly from §5$), 



dP{x,t) 
dt 



d_ 

dx 



D(x 



dP(x,t) 
dx 



with a spatially varying diffusion coefficient D(x) = k%T / (m~f(x)). 
dynamics of an ensemble of test particles is investigated on the basis of 
with their initial positions at x(t = 0) = 0. 



(6) 

The Brownian 
P), © and © 



2.2. Model viscosities 

Spatially varying viscosities may be realized with organic gradient materials, where for 
instance the composition or the degree of polymerization changes in space [33| [35] . In 
photorheological fluids, as another example, the viscosity can be tuned by a spatially 
varying illumination intensity [37] . A further class with the possibility of inhomogeneous 
viscosities are binary-fluid mixtures. Here spatial variations of the concentration of 
the two constituents may be driven by temperature modulations via the Soret effect 
[38] . If both constituents have sufficiently different viscosities the thermally induced 
concentration variations are accompanied by spatial viscosity changes. For materials 
with a strong Soret effect, quantified by the Soret coefficient St, small temperature 
gradients are sufficient to generate large concentration and therefore large viscosity 
gradients, so that the direct effects of temperature gradients can be neglected. Suchlike 
experimentally favorable large values of St can be achieved by shifting the mean 
temperature of the binary fluid in the one-phase region close to the critical temperature 
of the mixture, where a transition to the two-phase region takes place (HJ H2J H31 HI]. 

For the materials mentioned above, one can imagine a number of spatially varying 
viscosities. A simple example of a viscosity of experimental relevance, being asymmetric 
with respect to the particles' initial position at x = 0, is 

r](x) = 7] + A?7tanh(x/L) with r] > \Arj\ > , (7) 

which may be approximated for \x\ <^ L by 

r)(x) = r/o + fj\x with fji = Ari/L. (8) 

Other generic viscosities are symmetric with a minimum or maximum at the particles' 
starting point. As a representative model of these viscosities, showing a smooth change 
from the value t]q at the extremum to the bulk value rjoo, we choose 

Ar) 

V(x) = Voo+ : , | /T , K , (9) 
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with k > 0, the viscosity contrast A77 = f]o — Voo, and the characteristic length L. A77 < 
corresponds to a viscosity with a global minimum and An > to a maximum at x — 0. 

Equation ([9]) covers several limiting cases, which can be identified by expressing 
7]{x) in for |x| < L by a power series r/ = r/^ + Arf^2™ =0 (— l) n |f " around 
x = and for |x| > L by an asymptotic series 77 = 7700 + Ar/ X]^Lo(~ 1)™ III 
In distinct parameter ranges each series can reasonably be approximated by the leading 
contribution, which is feasible for the analytical considerations in section [3j 

For a pronounced minimum, i. e. r] <C t]^, two approximations for r)(x) are 
useful: ?7 ~ ?7o +r] 00 \x/L\ K for \x\ <C L respectively n ~ 7700 (1 — for 
|x| L. The viscosity reaches its constant bulk value rj^ for large values of \x\ and for 
Vo/Voo *C <C 1 the latter approximation evaluates to the power law 

r](x) oc \x\ K (10) 

that covers the x dependence of 77(2) rather well. The validity of this power-law can be 
extended over a wider range with decreasing ratio 770/7700- For a pronounced maximum, 
i. e. 7]q ^> f]oo, we approximate r](x) for \x\ <C L as 77 ~ 770 (1 — and for |x| ^ L 

as 77 ~ 7700 + 77o|x/L|~ K . Accordingly for 1 \x/L\ K <C 770/7700 the viscosity can again be 
simplified to a power law 

77(x) oc \x\~ K , (11) 

while for \x/L\ K ^> r/0/7700 a nearly constant viscosity with 77(2) « n^ results. The 
range, where the viscosity can reasonably be approximated by a power-law, increases 
with rising values of the ratio 770/7700 • 



3. Scaling arguments and analytical solutions 

For symmetric viscosity profiles, 77(|x|), we present a scaling analysis for the power law 
of the mean square displacement (x 2 ) oc t a and in limiting cases exact solutions of the 
Fokker-Planck equation ([6]). For asymmetric profiles we use a perturbation series to 
gain the time evolution of the first moments. 



3.1. Scaling arguments for symmetric viscosity profiles 

The mean square displacement of a particle (x 2 (t) ) is commonly used to characterize its 
random motion. For a constant viscosity, the analytical solution of equation (T5]) takes the 
well-known form [2]: (x 2 (t)} = 2kB f (e -7 * + 'yt — 1). Beyond a short period of ballistic 
motion in the regime of normal diffusion, t ^> 7 -1 , the mean square displacement 
increases linearly in time: 

{At)) = 2 lsl t = ^Li = Q i. (12, 

77?/y Sira 77 77 

In order to achieve progress by analytical calculations also for spatially varying and 
symmetric viscosity profiles 77 (x) = 77(|x|), we use (|T2l) with the replacement 77 — > 77(|x|) 
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and simultaneously the substitution |x| — > (x 2 (t)). The resulting expression, 

(At)} = Q f / 1 x , 

ri(^WW) 



is further analyzed for different regimes of the model viscosity (Q. 

Within a very short time interval, where (x 2 ) ~ holds, the model viscosity 
li^s/ ( x2 (^))) given by (Q can be approximated by t]q (regime I) and in the long time 
regime, i. e. for (x 2 ) — > oo, by 7]^ (regime III). In both cases one obtains with equation 
(TT5j) normal diffusion: 

( x 2 { t )) = -t and (x 2 (t)} = — t. (14) 

VO Voo 

Between these two regimes, one finds regime //with a power-law behavior rj (x) = fj 
respectively D(x) = D\x\~^ (see also ([6]), ( JTUl) and (lllj) ) for /3 = k > near a viscosity 
minimum and /3 = — n < near a maximum. This yields together with (fl~3j) 



(i2(i »= (3 ^4)^- (i5) 

which is the basis of the prediction of anomalous diffusion in regime II: 

(x2(t))= (ir * a * a= 2^- (i6) 

Near a viscosity minimum, (3 > 0, one has subdiffusion with a < 1 and near a viscosity 
maximum, /3 < 0, superdiffusion with a > 1. For j3 = —2 the scaling formula (1161) breaks 
down and the overdamped Langevin equation (jSJ) evaluates to x = x [2D + (2D) 1 ^ 2 £(£)] , 
which describes the so-called geometric Brownian motion. In this case the power law 
in equation ([T]) is replaced by an exponential time dependence of the second moment 
(x 2 (t)) oce 6 ^ @5] . 

During the crossovers between the three regimes the exponent a of the mean square 
displacement, (x 2 (t)) oc t a , varies as a function of time from 1 in regime I to 2/(2 + 0) 
in regime II and back to 1 in regime III. The location of regime II follows by translating 
the respective spatial ranges described in section [2T21 onto time via (fT3|) : 

2L 2 Vo (r) \ 2/K ^ f ^ (subdiffusion) ^ (17) 



Q \VooJ 2Q 
L 2 Vo <<; t ^ 2I^oc ( fjo\ V K (superdiffusion) . (18) 



2Q Q \Voo 

Therefore, the experimental accessibility of the anomalous regime II is enhanced, when 
the viscosity contrast \Ar]\ = \vo — Voo\ is enlarged. 

3.2. Fokker-Planck equation: Solutions for symmetric and asymmetric viscosity profiles 

For a power law approximation rj(x) oc \x\^ of symmetric viscosity profiles in regime 
II the diffusion coefficient is given by D(x) = D\x\~@. In this case, the Fokker-Planck 
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equation fl2]) with the initial condition P(x, 0) = 8(x) can be solved analytically for 
f3 > —2 by using the ansatz 

k ( lxl 2+/3 \ 

p ^ = wmj^{~jE(kw-) • (19) 

with k = 2±£ [D(2 + /3) 2 ]~V(2+/3) p-i(_L_) being a normalization constant. With this 
expression the second moment can be evaluated exactly to 

<x 2 > = \D(2 + /3) 2 tj 2/(2+/?) oc t 2 '^ , (20) 

which follows the identical power-law in time as in ( Fl6l) and supports the validity of the 
assumptions made for the scaling argument in section 13.11 The odd moments vanish 
because of the ±x symmetry and the even moments deviate for (3^0 from the case of 
normal diffusion. As an example we consider the kurtosis K, which is a quantity for the 
distribution's peakedness and defined as 4th central moment divided by the 2nd central 
moment [8j. For ( [19]) it takes the form K = T^^aj -^(2+^) ^~ 2 (l+^)' indicating a 
normal distribution (K = 3, i. e. mesokurtic) for (3 = 0, a distribution with a high 
sharp peak (K > 3, i. e. leptokurtic) for — 2 < j3 < 0, and a flat-topped distribution 
(K < 3, i. e. platykurtic) for j3 > 0. 

To gain insight in the diffusion for a asymmetric linear viscosity profile as in 
equation (|E]), respectively for the corresponding diffusion coefficient, D(x) = D(r] + 
?7ix) _1 with D = we derive the equation for the n-th moment, 

d ( x n ) = n(-^-D(x)x n - 1 \ , (21) 



dt \dx 

as a function of D(x) by a double integration by parts of the Fokker-Planck equation ([6]). 
In order to proceed we expand on the right hand side of equation (12 lj) the term inside 
the brackets with respect to powers of x. The resulting system of coupled differential 
equations for the moments (x n ) is solved in the limit fjix/riQ C 1 by a power series 
ansatz in fji. This leads together with the initial conditions (x n )(0) = to the following 
approximations for the first three moments: 

(x) = -^t-^!t 2 + 0(^), (22a) 
Vo Vo 

(x 2 )= —t + -±^t 2 + '\ t 3 + 0( V f), 226 
Vo Vo 3 Vo 

, a . 12771.D 2 2 UOffiD 3 o m . , N . . 

(x 3 ) = '\ — t 2 §— t 3 + 0{r,l) . (22c) 

Vo Vo 

According to ( 122 &|) anomalous diffusion can be expected to occur for linearly varying 
viscosity profiles, too. Further, the onset of anomalous diffusion is governed by the ratio 
771/7/0: The larger this ratio the sooner the transition to anomalous diffusion takes place. 
Asymmetric viscosity profiles lead to a mean particle drift towards the region with lower 
viscosity as indicated by (|22aJ). It is also a remarkable property of each of the three 
moments that the sign of the leading contributions does not alternate. 
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4. Numerical results 

The results of the previous section are substantiated in the following by comparisons with 
numerical simulations of the Langevin equation (j2J) and the Fokker-Planck equation ([6]) 
for different spatially varying viscosities. 



0.4 




Figure 1. The probability density function P{x) at t = 10 is given for -q(x) = ?7o +771^1 
with r/o = 10 and fj = 10 in (a) and fj = 100 in (b). The solid line (FPE) is obtained as 
a numerical solution of (|6|) with a narrow initial distribution and the step function (LE) 
by integrating ([2]) for 10 6 independent trajectories with the starting point at x = 0. 

The Fokker-Planck equation ([6]) with a spatially varying damping is derived from 
(J2J) under the assumption t ^> 7 -1 [M]. For the symmetric viscosity 77(2) = rj + fj\x\ 
with the minimal value r] = 10 and two different slopes fj = 10 and fj = 100, the 
distribution P(x) is determined in figure [Hat the time t — 10 3> 7~ iT1 = 0.1 by solving 
either the Fokker-Planck equation fl6]) (black solid curve) or by integrating the Langevin 
equation ([2]) for an ensemble of 10 6 particles (red step function). The inequality t 3> 7 _1 
is not always fulfilled as good as in figure [U but the solutions of ([6]) are also for shorter 
times and stronger variations of r](x) often in good agreement with the solutions of (J2]). 

There is a further aspect to be drawn from the profiles of the probability distribution 
P(x) in the anomalous regime. After the short normal diffusion in regime I, P(x) at 
t — 1 deviates in regime II from a Gaussian profile as shown in figure EJ In part (a) 
the red curve of P(x) is determined numerically by solving fl6]) for the spatially varying 
viscosity r)(x) = t] + x 2 with a small minimal viscosity 770 = 0.01. It is compared 
with the analytical solution of the Fokker-Planck equation in ( f!9l) for rj = x 2 (blue 
dashed line) and we find a nearly perfect agreement between both approaches in figure |2j 
Deviations grow by increasing the viscosity minimum 770 in our simulations. Both curves 
in figure [2]^a) deviate, however, significantly from a Gaussian distribution (green dotted 
line) with the same second moment. As predicted in section 13.21 the distribution is 
platykurtic for > 0, which can be explained as follows. Due to the small minimal 
viscosity 770 at the starting point at x = 0, particles quickly diffuse away. This reduces 
the probability to find a particle at x = and simultaneously enhances P(x) in a wider 



Anomalous diffusion in viscosity landscapes 



8 



neighborhood of the minimum of 7(2;), compared to the case of a constant damping 
with a Gaussian distribution. On the other hand, a viscosity increasing strongly as a 
function of |x| impedes quick particle diffusion away from the minimum. This reduces 
the probability to find a particle at larger distances compared to a Gaussian distribution. 



0.4 




Figure 2. The probability distribution P(x) at t = 1 for two different viscosities: In 
(a) T)(x) = 0.01 + x 2 with a minimum at x = and in (b) rj(x) = (0.001 + | £c | ) 1 with a 
maximum at x = 0. The red line in each part is obtained by numerical integration of 
(|6]) and the nearly identical blue dashed curve is the analytical solution given by (TH)]) 
for r](x) oc x 2 in (a) and rj{x) cx | a; | 1 in (b). In both parts P{x) is compared with a 
Gaussian profile (green dotted line) having the same second moment. 

With a viscosity maximum at the starting point of the particles the situation 
is similar, as illustrated by the distribution in figure |2|b). Here the numerical 
solution P(x) of the Fokker-Planck equation (red curve) is obtained for the viscosity 
i](x) = (0.001 + I I ) 1 and is compared with the analytical solution given by ( TT9|) for 
t](x) oc (blue dashed curve). Only in the vicinity of the viscosity maximum at 

x = both approaches differ slightly. Again we find clear deviations from a Gaussian 
distribution with an identical second moment as described by the green dotted curve in 
figure |2^b). Since the Brownian dynamics is significantly reduced near the maximum of 
the damping at x — 0, the particles move much slower away from the peak. This causes 
an enhancement of P(x) around the maximum, compared to a Gaussian distribution. 
Thus the distribution is leptokurtic, in accordance with our predictions in section 13.21 

For a further analysis of the Brownian motion the mean square displacement 
(x 2 (t)) of the particle positions is calculated by integrating the Langevin equation (T2]) 
numerically for an ensemble of particles starting at x — 0. Figure |3] exemplarily shows 
(x 2 (t)) as a function of time for two different spatially varying viscosities. In each case 
and independent from the particular viscosity landscape, the mean square displacement 
is characterized by the initial regime I of normal diffusion, marked by the dotted lines. 
Beyond this regime I, anomalous diffusion occurs with exponents a 7^ 1, as long as 
the viscosity, experienced by the particles, does not reach a non-vanishing constant 
plateau value, as for instance for large \x\ in (JSJ). In the vicinity of a viscosity minimum, 




Figure 3. The mean square displacement (x 2 (t)) as obtained by a numerical 
integration of © for two different model viscosities r](x). In (a) 77(21) = 10 + |x| 
has a minimum at x = 0, whereas in (b) one has for |x| < 1 the plateau tj(x) = 20, 
which is continued for \x\ > 1 by the decaying function r/(x) = 20|x| -3 / 4 . The dashed 
and dotted lines in both parts indicate linear fits to the numerical data, whereby the 
slope a is explicitly shown. 



subdiffusive Brownian motion with exponents a < 1 is found, as exemplified in figure 
[3(a) for T](x) = f]o + fj\x\ with rj = 10 and fj = 1. At long times the mean square 
displacement (a; 2 ) scales in the range of the dashed line with t a and a = 0.69, in good 
agreement with the scaling result a = 2/3 predicted by formula Q16j) with /3 = 1. The 
viscosity used in figure [3(b) has in a small range |x| < la constant maximum, rj(x) = 20, 
and then decays in the range \x\ > 1 according to the power law rj{x) = 20|x|~ 3 / 4 . A fit 
to the numerical data in the limit of long times yields in the range of the dashed line 
in figure [3(b) the exponent a = 1.58, deviating only slightly from the scaling prediction 
a = 1.6 obtained from f fl6|) for /3 = —0.75. 

The deviations between the exponent obtained in simulations and the exponent 
obtained via scaling arguments decrease at sufficiently long times either by reducing 
the minimal viscosity ?]o as in figure [3(a) or by shortening the plateau in figure [3(b). 
The scaling regimes I (normal diffusion) and II (anomalous diffusion) are separated 
by a transition period, where the crossover from a ~ 1 to a ^ 1 takes place. This 
transition period may extend over several decades in time as for instance in figure [3(a) 
with 77/770 = 0.1. Together with an exponent (3 = 1 this causes a transition period 
lasting over more than three decades. In contrast, steep, nonlinear viscosity gradients 
obtained for (3 > 1 and (3 < delimit the intermediate regime [see figure [3(b)] . 

For a further examination of the scaling prediction in the anomalous regime II, 
different spatially varying viscosities were investigated numerically and in terms of the 
analytical scaling argument given by (TIB"]) . As can be seen in table 1, the numerically and 
analytically obtained exponents a are in excellent agreement for several viscosities. In 
general, in the neighborhood of a viscosity minimum one obtains subdiffusive behavior, 
whereas in the vicinity of a viscosity maximum superdiffusive behavior is observed. 



Anomalous diffusion in viscosity landscapes 



10 



viscosity rj(x) 




a (scaling) 


a (numerically) 


n m x M 4 

U.Ul \X\ 


A 
-± 


1 l'\ 

1/6 


u.ooo 


n ni i 1 ™l 2 
U.Ul + \X\ 





1 /o 
1/2 


U.OUl 


0.01+ a; 


1 


2/3 


0.666 


0.01 + \x\^ 2 


1/2 


4/5 


0.799 


0.01 + l^l 1 / 4 


1/4 


8/9 


0.886 


(0.01 + |^| 1 /2)-i 


-1/2 


4/3 


1.330 


(0.01 + i^i 2 / 3 )- 1 


-2/3 


3/2 


1.494 


(0.01 + Ixl)- 1 


-1 


2 


1.976 


{0.01 + \x\ 4 / 3 )- 1 


-4/3 


3 


2.963 



Table 1. Chart of anomalous exponents 
a for spatially varying viscosities rji. x ) 
obtained by the scaling formula (fl~6|) and 
via numerical solutions of ([2]) by fitting 
the resulting mean square displacement 
in regime II by the power law cx t a . We 
find a good agreement in both generic 
cases: For subdiffusion (a < 1) with a 
minimum of rj{x) and for superdiffusion 
(a > 1) with a maximum of r}{x). 



For the viscosities in table 1 with a small value rjo at the minimum or a large value 
at the viscosity maximum one finds according to ([17]) and ( TT8l) an early onset of the 
intermediate anomalous regime II. We would like to emphasize, that the viscosities 
discussed up to now are approximations of (jHJ), which we used to demonstrate the 
good agreement between our scaling results and numerical simulations in the anomalous 
regime II. For the viscosity ([9]) regime II is also limited in time from above as described 
in the previous section by the characteristic times given in (II 7p and ([IF 



0.9 



0.8 



0.7 



(a) 



° Vo = 1 
x Vo = 10 
+ 770 = 100 



-1 1 
logr) 




-2-10 1 
log(?7/?7o /2 ) 



Figure 4. The scaling exponent a of the mean square displacement, (a; 2 ) cx t a , is 
plotted in (a) for the symmetric viscosity rj(x) = t/q + }jj|a;| as a function of the slope fj 
for three different values of 770 = 1, 10, 100 and in part (b) as a function of the ratio 
v/Vq^ 2 - a i s obtained by fitting (a; 2 ) cx t a over a fixed time interval < t < 10 2 to 
numerical data created by integration of equation 



To investigate the crossover behavior between regime I and II, the Langevin and the 
Fokker-Planck equation were solved for the viscosity 77 = 770 + fj\x\ up to the fixed time 
T = 100 and for different values of r] Q and fj. With this approach we mimic the limited 
time range available in experiments for detecting the anomalous regime. The numerical 
curves of the mean square displacement (x 2 (t)) are fitted over the whole range < t < T 
by the power law t a and the resulting exponent a is shown in figure H] as function of the 
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viscosity. The deviation of a from a = 1 grows, when the anomalous diffusion regime 
II occupies an increasing part of the interval [0,T], i. e. when the inequality fj\x\ > t] Q 
holds for an increasing part of [0, T\. The two terms f}\x\ and rjo in the denominator of 
( TT3|) become nearly equal at the time t t ~ Vo/v 2 - Hence, for decreasing values of t t the 
range of anomalous diffusion in [0, T] increases as illustrated by the trend in figure H^a). 
On the other hand, the anomalous fraction and therefore a can be kept constant by 
keeping r\ 2 jr)\ constant, which is shown by figure HJ^b). 

The behavior of a in figure H] can also be derived via ( fl3l) by using 77 = r]Q + f](x 2 ) 1 ^ 2 . 
The resulting nonlinear equation for the second moment, (x 2 )(r] + ^(x 2 ) 1 / 2 ) = Qt, can 
easily be solved. From this solution the exponent a(t), which describes the local slope 
along curves as in figure [31 can be calculated by a(t) = d\og((x 2 )) / d\og(t) . Its average 
a = Jq dt a(t) corresponds to the dashed lines in figure HJ which agree surprisingly 
well with the full numerical results obtained by ([6]). We find visible deviations from 
the full solution only if the viscosity shows very small values close to zero either at the 
minimum or for \x\ 3> L, which, however, is unlikely in experiments. Therefore, the 
exponent obtained via the scaling relation f fl3|) may be useful in analyzing and fitting 
experimental results. 

A further experimental relevant model viscosity is described by the asymmetric 
function ([7]), which may be approximated in the range \x/L\ <C 1 by the linear 
dependence in equation (IE]). The essential qualitative difference to the previous model 
viscosities is its asymmetry with respect to the starting point of the particles at x — 0, 
which causes non-vanishing moments (x) and (x 3 ) as given in (122 al) and (122 di) for short 
time intervals. The quality of the approximate solutions in ( I22q|) - (l22cl) is investigated 
in figure [U where in part (a) the leading contribution to the drift of the mean value in 
(I22ap (dashed lines) is compared to the first moment of the numerical solution of ([6]) 
(solid lines). Despite the finding that the differences between both approaches increase 
with fjit, the results are still in good agreement since the deviations do not exceed 0.25% 
for fji = 0.25 and 4.5% for fji = 1, even at t = 10. 

The diffusion in an asymmetric viscosity profile (IE]) becomes anomalous too, as 
indicated by the contributions oc t 2 and oc t 3 to the approximate second moment 
in equation ( 12251) . Again, the deviations are small between the full numerical result 
of (x 2 (t)) [solid line in figure EJ^b)] and the approximation ( 12261) [dashed line in 
figure 0(b)] . We point out, that anomalous diffusion still persists, when the drift of 
(x) is subtracted from the particle dynamics as illustrated by the following formula: 
(x 2 ) - (x) 2 = 22 t + ^ t 2 + 0(fjf). With the help of the third moment (12211 one can 
also define the skewness 5* of a probability distribution [8] , which is given in our case at 
the leading order by S = —(3f]i\/Dt)/ Since S is negative for fji > 0, the left 
tail of the probability distribution becomes more important. 

All three aspects, the non-vanishing drift, the skewness of the particle distribution 
and the anomalous diffusion caused by an asymmetric monotonous viscosity profile are 
recaptured in numerical solutions of the Langevin equation ([2]) and the Fokker-Planck 
equation (J6j) by integrating both for the viscosity ([7]) up to time scales beyond the 
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Figure 5. The first two moments of the numerical solutions of © (solid lines), (x) 
in part (a) and (a; 2 ) in (b), are compared with the leading order contributions to the 
approximations given by (|22q|) respectively (|226[) (dashed lines) for a linear viscosity 
profile 97(0;) = 7/0 + fjix with r/o = 10 and different values of 771 = 0.25, 0.5, 1. 
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Figure 6. The probability distribution P(x) at t = 1 in (a) and the mean square 
displacement (x 2 ) in (b) are obtained by integrating for the asymmetric viscosity 
([7]) with ?7o = 10 and different values of A77 = 0, 5, 9. The inset in (b) shows the 
three regimes of (x 2 (t)) for 770 = 10 and A77 = 9. In the intermediate regime II at the 
largest slope (dashed line) the scaling exponent a is determined by the fit (x 2 ) cx t a to 
the numerical data. In (b) the resulting a in regime II is given as function of A77/770: 
The red solid line belongs to the evolution of (a; 2 ) and the dashed line to the variance 



,2\ 



validity range of equations ( I22ap -( )22cj) . The asymmetry of the resulting probability 
distribution is shown in figure [6]^a) for increasing values of An, where the left tail is 
the more important one as predicted by the skewness S of the distribution introduced 
above. The time dependence of the mean square displacement (x 2 (t)) is shown for the 
viscosity ([7]) by the inset in figure Mb) for L = 1, rjo = 10 and Arj = 9 and illustrates 
all three temporal diffusion regimes. During the early temporal regime I with \x\ L 
and in the long time limit with |x| ^> L (regime III) one has normal diffusion with the 
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exponent a = 1 as indicated by the dotted lines in the inset of[^b). In contrast, during 
the intermediate regime II with |x| ~ L one has superdiffusion with a = 1.543 > 1 as 
indicated by the dashed line. If the numerical data of the mean square displacement, 
obtained for t] = 10 and Arj = 9, are fitted by (x 2 ) oc t a in the range of its largest slope 
in regime II we find that a is independent of the length L as indicated by the values 
a = 1.538 (L = 0.2), a = 1.543 (L = 1), and a = 1.544 (L = 5). However, the exponent 
a varies strongly as a function A 77/770 in figure Mjo). As predicted by the analytical 
expressions m (T22fflj) and ( 1221 the mean square displacement (x 2 ) [see also figure EJ^b)] 
and the variance (x 2 ) — (x) 2 show anomalous diffusion behavior as well. The solid line 
in figure EJ^b) describes the exponent a of the mean square displacement and the lower 
lying dashed line is the exponent related to the variance (x 2 ) — (x) 2 . For both quantities 
the exponent a may take even larger values in the limit Arj/rjo — > 1. The mean value 
(x) drifts also for the viscosity profile (I7j). For small values of t and x we find a perfect 
agreement with the linear time dependence given by the leading contribution in (122 fej) . 
but its dependence becomes more complex in the range of larger values of A77 and in 
the long-time limit. 



5. Conclusions 



We have identified three different diffusion regimes for mesoscopic Brownian particles 
in spatially varying viscosities, which were analyzed by three approaches: Firstly, by 
scaling arguments applied to the expression of the mean square displacement, secondly 
by simulations of the corresponding nonlinear Langevin equation (T5]), and thirdly by 
solving the related Fokker-Planck equation fl6]) either numerically or, in limiting cases, 
even analytically. For an ensemble of particles starting at a viscosity extremum a short 
regime of normal diffusion is found where the mean square displacement (x 2 ) oc t a 
scales with the exponent a = 1. Beyond this regime Brownian particles experience 
considerable changes of the viscosity along their trajectories, which leads to anomalous 
diffusive motion in an intermediate temporal regime. Near a minimum of the viscosity 
the particle dynamics becomes subdiffusive with an exponent a < 1, whereas it becomes 
super diffusive with a > 1 in the vicinity of viscosity maxima. In the long time limit, 
when the Brownian particles explore the whole range of a viscosity variation, as for 
instance described by or for a periodically varying viscosity, the particles experience 
a mean viscosity and therefore normal diffusion is found. 

In the intermediate anomalous diffusive regime the particles' probability 
distribution is non-Gaussian. In the case of a viscosity minimum the particle distribution 
function is - compared to a Gaussian with the same second moment - reduced at its 
maximum, increased at intermediate distances from its maximum, and reduced again at 
large distances. The opposite happens near a maximum of the viscosity: The particle 
distribution function is enhanced at the starting point of the particles and reduced at 
intermediate distances from the initial position. For an ensemble of particles starting in 
the range of a linearly varying viscosity profile, which is asymmetric with respect to the 
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starting point, one finds an asymmetric particle distribution reflecting itself in a non- 
vanishing skewness. In addition, the first moment of the distribution drifts as function 
of time into the direction of decreasing viscosity and the particles show superdiffusive 
behavior. However, superdiffusivity is only partly caused by the particle drift: The 
exponent of the superdiffusive motion reduces only slightly towards 1 when the mean 
drift is subtracted from the particle motion but still shows clearly superdiffusive behavior 
with a > 1. 

The subdiffusive behavior in an intermediate regime in the neighborhood of a 
minimum of the particle damping shares similarities with the Brownian dynamics of 
a single segment of a flexible polymer, the so-called Rouse dynamics. Also in this case 
one has normal diffusion of the segments on a short as well as on a long time scale. In 
the intermediate range, where the mean square displacement of the segment is of the 
order of the polymer-coil diameter, one finds the subdiffusive behavior (x 2 ) oc t 1 / 2 [32J. 

In contrast to most of the well-known examples showing anomalous diffusion, one 
can imagine for systems suggested in this work both types of anomalous diffusion, sub- 
and superdiffusion. In photorheological materials [37], for example, either a maximum 
or a minimum of the viscosity as well as monotonously varying viscosities can be induced 
by an appropriate spatial variation of the illumination strength. In some binary mixtures 
the sign of the Soret effect changes as a function of the mean temperature of the mixture 
[4"6] . The two possible signs may be used to attract either the lower viscous component 
to the locally heated area or the higher viscous component. In glass forming polymer- 
mixtures the relation between the local composition and the viscosity can be a nonlinear 
function [39l HQ]. Materials with strong variations of the local viscosity are especially 
favorable to observe anomalous diffusion, as discussed in this work. 

Another interesting subject are heated particles in a binary fluid mixture, 
undergoing Brownian motion. An example are light absorbing particles in a transparent 
binary fluid-mixture with an upper miscibility gap such as for instance n-Butoxyethanol 
and water [H], where the temperature may be driven close to the transition temperature. 
The inhomogeneous temperature field around a particle sets in much faster than related 
changes of the concentration and the viscosity in the neighborhood of the particle, such 
that the Brownian particle experiences temporally its own induced viscosity variations. 
The delayed dynamics between temperature and viscosity field may lead to interesting 
memory effects that are also prone to cause anomalous diffusion and will be investigated 
in forthcoming work. 
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